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ABSTRACT 



We report precise Doppler measurements of GJ317 (M3.5V) that reveal the 
presence of a planet with a minimum mass Mp sini =1.2 Mjup in an eccentric, 
^ ! 692.9 day orbit. GJ317 is only the third M dwarf with a Doppler-detected Jo- 

vian planet. The residuals to a single-Keplerian fit show evidence of a possible 
second orbital companion. The inclusion of an additional Jupiter-mass planet 
^ ; (P « 2700 days, Mp sini = 0.83 Mj^p) decreases ^/xl from 1.42 to 1.11, and 

reduces the rms from 12.5 m s~^ to 6.32 m s~^. A false-alarm test yields a 
1.1% probability that the curvature in the residuals of the single-planet fit is due 
^ ! to random fluctuations, lending additional credibility to the two-planet model. 

O '■ However, our data only marginally constrain a two-planet fit and further moni- 

^ i toring is necessary to fully characterize the properties of the second companion. 

I To study the effect of stellar mass on giant planet occurrence, we measure the 

O I fraction of stars with planets in three mass bins comprised of our samples of 

Q I M Dwarfs, Solar-mass stars, and intermediate-mass subgiants. We find a pos- 

itive correlation between stellar mass and the occurrence rate of Jovian planets 
I within 2.5 AU. Low-mass K and M stars have a 1.8 ± 1.0% planet occurrence 

^ I rate compared to 4.2 ± 0.7% for Solar-mass stars and 8.9 ± 2.9% for the higher- 

mass subgiants. This result indicates that the former F- and A-type stars with 
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M* > 1.3 M0 in our sample are nearly 5 times more likely than the M dwarfs to 
harbor a giant planet. Our analysis shows that the correlation between Jovian 
planet occurrence and stellar mass remains even after accounting for the effects 
of stellar metallicity. 

Subject headings: techniques: radial velocities — planetary systems: formation — 
stars: individual (GJ317) 



Introduction 



A planet host star inherits its fundamental characteristics from the same disk material 
that forms its planets. Studying the relationships between the observed occurrence rate 
of giant planets as a function of the characteristics of their host stars therefore provides 
crucial tests of planet formation theories. The two most fundamental properties of stars are 
mass and chemical composition. A number of detailed spectroscopic analyses of nearby stars 
have revealed a strong correlati on between the metallicity of star s and the likelihood tha t 
they harbor detectable planets (iGonzalezI 119971 : ISantos et al.l 12004 iFischer fc Valentill2005l ). 
This finding can be understood in the context of the core accretion model: increasing the 
metallicity of star/disk system increases the surface density of solid particulate matter, which 



leads to an enhanced growth rate for protoplanetary cores ( jida fc LinI l2004l : iKornet et al. 
2005l : llda fc Lh]l2005ah . 



Another way to enhance the surface density of solid mate rial in the midpla ne of the 
protoplanetary disk is to increase the total mass of the disk ( llda fc LinI l2005bl ). If disk 
masses scale with the mass of their central stars, then there should be an observed correlation 
between planet occurrence and stellar mass. The relationship between stellar mass and planet 
formation rates has b een explored theoretically in the context of the core accretion model. 
Laughlin et al.l (120041 ) simulated the formation of planets in disks around low-mass stars 
and found that the lower surface densities and longer orbital time scales in the disks around 
M dwarfs impede the growth of Jupiter-mass planets. By the time the critical core mass is 
reached (~ 10 Me), the s upply of disk ga s is exh austed due to accretion onto the central star 
and photo-evaporation. iLaughlin et al.l (120041 ) therefore predict an abundance of "failed" 
gas giants with masses comparable to Neptune and a much lower frequency of Jupiter-mass 
planets. 

The mass of the c entral star also influences th e radial extent of the region in which proto- 
planetary cores form. [Kennedy fc Kenyoru (120071 ) account for the evolving luminosity of the 
central star to model the influence of stellar mass on the location and size of the core-forming 
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region in circumstellar disks. Similar to iLaughlin et al.l (120041 ) a ndllda &: LinI (2005bl ) , they 
predict that M dwarfs should exhibit a deficit of giant planets. Kennedy fc KenyonI (120071 ) 
also predict that the fraction of stars with Jovian planets should increase with stellar mass 
up to a peak of ^ 20% for = 2.5 Mq. 

The prediction of a lower frequency of Jupiter-mass planets around low-mass M dwarfs 
is in agreement with the current observational data. About 300 nearby M dwarfs are cur- 
rently monitored by various Dqppler surveys, and only 5 have been discovered to harbor one 
or more planets (IDelfosse et al 



2005b: Butler et al 



2006 
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19981: iMarcy et al.lllQQsl. I2OO1I : ISndl et al.l[2003l : ISonfils et al" 



Bonfils et al 



20071 ). Three additional planets have been discov- 



ered orbiting distant (~ 1 kpc) low-mass sta. r s by gravitational le nsing surveys (IBond et al. 
20041 : IGouM et ahlbood : iBeaulieu et allbood : iBennett et aUbooeh . 



The majority of the Doppler-detected planets around M dwarfs are significantly smaller 
than Jupiter, with minimum masses {Mp sin i ) comparable to the masses of Neptune and 
Uranus (5-23 M®). The sampl e of low-mass, Doppler-detected planet s consists of the 22.6 ± 



1.9 M® planet around GJ436 (iButler et al.ll2004J: iManess et al.ll2007t l: the inner, Mp sin i 



5.9 planet around GJ 876 (IRivera et al.l l2005l ) : the Mp sin i = 11 planet orbiting 
GJ674 (IBonfils et al.ll2007l ): and the triple syste m orbiting GJ 581, which consists of p lanets 
with minimum masses 16.6 M®, 5 and 8 ( IBonfils et al.ll2005bl : lUdry et al.ll2007l ). The 
kinship between these extrasolar "super Earths" and the ice giants in our Solar System was 
confirmed wi th the recent disco very that the Neptune-mass planet around GJ436 transits its 
central star (iGillon et al.l 120071 ). The transit light curve provides an absolute measurement 
of the planet's mass and radius, suggesting an internal structure very much like Neptune, 
with a rocky core likely surrounded by a thick water layer and thin gaseous envelope. 

To date, there are only two nearby M dwarfs kno wn to harbor Jupiter-mass compa nions: 
GJ876 jMarcv et al.lll998l . I2OO1I : iRivera et al.ll2005h and GJ849 (|Butler et al.ll2006ah . And 
so far not a single "hot Jupiter" (P < 10 days) has been discovered around an M star, in 
contrast to the 1.2 ± 0.8% occurrence rate of hot Jupiters around Sun-like stars (lEndl et al. 



2OO3I : iButler et al.ll2006bl ). The relatively small number of Jovian planets around M dwarfs 
is not due to a decreased sensitivity of Doppler techniques for these stars. For a planet 
of a given mass and orbital period, the amplitude of a host star's refiex velo city scales as 



2/3 

K PC M ^ , making planets easier to detect around stars with lower masses. IButler et al. 



( l2006al ) showed that the frequency of giant planets is 2-3 times higher arn ong Solar-mass 
stars compared to M dwarfs (See also iLaws et al.l l2003l : IBonfils et al.l 120071 ). However, the 
uncertainty in the estimated planet occurrence rate for M dwarfs is large due to the small 
number of target stars and planet detections. 



Another obstacle that has so far hindered a study of the effects of stellar mass on planet 
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formation is the limited range of masses spanned by Doppler-based planet searches. For 
example, the Solar-mass, FGK dwarfs that comprise the bulk of the California and Carnegie 



Planet Search (CC PS) span stellar masses from 0.8 Mq to 1.2 Mq (jValenti &: Fischerl 12005 



Takeda et al.l 120071 ) . Much more leverage can be gained by measuring the planet occur- 
rence rate around higher mass F- and A-type stars (1.3 < M* < 3.0 Mq). Unfortu- 
nately, intermediate-mass, main-sequence stars are poor precision Doppler targets. Stars 
with spectral types earlier than F8 tend to have ro tationally broadened absorption fea- 
tures (Ido Nascimento et al.ll2003l : iGalland et al.l 120051 ) . have fewer spectral lines due to high 
surface temperatures, and display a lar ge amount of excess velocity s catter due to surface 
inhomogeneities and pulsation ("jitter": ISaar et al.lll998l : IWrightll2005l ). 



We have addressed the limited range of stellar masses in Doppler-based planet searches 
by conducting a search for planets around stars with masses bracketing those of the Sun-like 
stars in the main CCPS sample. At higher stellar masses we are conducti ng a search for 



planet s arou nd intermediate-mass subgiants at Lick and Keck Observatories (jjohnson et al. 



2006al . 120071 ). Subgiants are evolved stars that have cooler surface temperatures and lower 
rotational velocities than their A and F type main-sequence progenitors (Tg// ~ 5000 K 
compared to T^jf > 6000 K, and K-otsini < 5 km s~^ versus K-otsini ^ 50 km s~^, respec- 
tively). The spectra of subgiants therefore have an increased number of narrow absorption 
lines required for precision Doppler measurements, making them ideal proxies for A and F 
type main-sequence stars. 

At the low-mass end we have been monitoring a sample of 147 late- K through M 



dwarfs (M, < 0.6 M^) as part of the NASA Keck M Dwarf Planet Survey (IButler et al. 



2006bl : iRauscher &: Marcyll2006l ). The 7-year baseline of our M dwarfs survey, together with 
the 3-year duration of our subgiants planet search now provides an excellent opportunity 
to measure the relationship between stellar mass and planet occurrence for wide range of 
stellar masses and comparable detection characteristics. Here, we report the detection of a 
Jupiter-mass planet in a 1.897 yr orbit around the M3.5 dwarf GJ 317. We present the stellar 
characteristics of the host star in § [2l In § [3] we present our observations and orbit solution, 
including an assessment of a possible second Jovian companion in the system. In § |4] we 
incorporate this latest M dwarf planet detection into a detailed analysis of the relationship 
between stellar mass and the occurrence rate of giant planets. We conclude in § [S] with a 
summary and a brief discussion of our results. 
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2. Properties of GJ 317 



GJ317 (LHS2037, L 675-081), is a nearby M3.5 dwarf among the 147 low-mass stars 
in the NASA Keck M Dwarf Survey. The Simbad database lists an apparent magnitude 

V = 13.0. However, this magnitude is ~ 1 mag faint er than the values listed in several other 
catalogs: The Gliese-Jahreiss Catalog (GJ Catalog: iGliese fc Jahreifilll99ll ) gives V = 12.0, 
the Carlsberg Meridian Catalo gs list V = 12.03, and the Catalogue of Stellar Spectral 
Classifications lists V = 11.98 (jSkiflg |200J) . We adopt the mean of these measurements, 

V = 12.0. The GJ Catalog also lists B — V = 1.52 and a trigonometric pa rallax tt = 



101.3 ± 27 mas. The GJ Catalog color agrees well with the color measured by iReid et al. 
(B — V = 1.5 3; I2OO2I ). and the listed parallax is consistent with the value measured by 
WooUeyl (Il970l ) (tt = 116 ± 12 mas). We adopt the mean of the two trigonometric parallax 
measurements, tt = 109. ± 20. mas. This value is 1-a larger (based on the quadrature- 
sum of the uncertainties) than the "resulting parallax" vr = 87± 17 mas listed in the the GJ 
Catalog, which is estimated from the star's broadband colors and spectral type. Our adopted 
trigonometric parallax yields a distance of 9.17 ± 1.7 pc and an absolute visual magnitude 
Mv = 12.2. 

For its B — V color GJ317 is unusually faint , lying f» 1.8 mag below the mean Hip- 
parcos main sequence as defined by IWrightl (120041 ). Its location in the H-R diagram sug- 
gests that GJ317 is an extremely metal-poor subdwarf. However, the star's 2MASS in- 
frared magnitudes together with the K-band photometric metallicity-luminosity calibra- 
tion of iBonfils et al.l (l2005a( ) suggest that GJ317 has a metallicity consistent with solar, 
[Fe/H]= —0.23 ± 0.2. This metallicity is much higher than the value expected given the 
star's position in the H-R diagram. 

We compared our iodine-free template spectrum of GJ317 to a template s pectrum of a 



similar M dwarf , GJ849 (M3.5, V = 10.4, B-V = 1.52, [Fe/H]= +0.16 ± 0.2. iButler et al. 
(j2006al )). Figure [1] shows a comparison of a small portion of the template spectra, which 
shows that these two stars are remarkably similar. Since the line depths are nearly identical 
we can safely conclude that GJ317 is not a metal-deficient subdwarf. As an additional 
check, we compared the 2MASS J — K and H — K colors of GJ317 to those of GJ849. 
The infrared colors of the two stars agree to within 0.001 and 0.006 mag, respectively, again 
suggesting very similar spectral characteristics. 

We can imagine several possible explanations for the abnormal location of GJ317 in 
the H-R diagram. The first possibility is that the parallaxes listed in the literature are 
systematically incorrect, and GJ317 is actually significantly further away than the reported 
parallaxes suggest. Similarly, the star's reported apparent V magniude may be too faint. 
Another possibility is that the star's B — V color is too blue by ~ 0.1 mag. Finally, the 
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Fig. 1. — Comparison of the spectra of GJ317 (black line) and GJ849 (gray line), which has a similar 
B — V color, spectral type and metallicity ([Fe/H]= +0.16 ± 0.2). The spectrum of GJ849 has been shifted 
up by 0.15 for clarity. Since the spectral lines of the two stars have very similar depths, we conclude that 
GJ317 cannot have a metallicity much different than that of GJ849. GJ317 is therefore not a metal-poor 
subdwarf, despite its apparent location 2 mag below the mean Hipparcos main sequence. 



star may be obscured by 1.5 mag of gray extinction. We are planning follow-up photometric 
monitoring during the next observing season to further investigate the latter three scenarios. 

We used the 2MASS infrared photo r netry of GJ317, together with the K-band mass- 
luminosity calibration of iDelfosse et al.l (120001 ) to estimate a stellar mass M^, = 0.26 ± 
0.04 Mq. This mass estimate agrees w ell with the . 22 ± 0. 05 Mq value predicted by the 
V-band mass-luminosity relationship of iBonfils et al.l (l2005al ). We adopt the mean of these 
two estimates, = 0.24 Mq. The properties of GJ317 are summarized in Table [2l 
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Fig. 2. — Radial velocity time series for four stable M dwarfs in our Keck Doppler survey demonstrating 
a Doppler precision of 3-4 m s~^ over the past 8-9 years. 



3. Observations and Orbital Solutions 

We have been monitoring GJ317 with the Keck I 10 m telescope for 7.4 years. We 
obtained high-resolution spectra using the HIRES echelle spectrometer (IVogt et al.lll994l ) 
with an iodine cell mounted directly in front of the entrance slit. The dense set of molecu- 
lar absorption lines provide a robust wavelengt h fiducial, as well as information about th e 
shape of the spectrometer instrumental profile (IMarcy fc Butlerlll992l : IValenti et al.lll995l ). 
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The Doppler shift is m easured from each star-plus-iodine observation using the modehng 
procedure described by iButler et al.l (119961 ) . Figure [2] shows velocity measurements of four 
stable M Dwarfs, demonstrating our long-term Doppler precision of 3-4 m s~^. 

A total of 18 Doppler measurements of GJ317 are listed in Tabled] and shown in Fig- 
ure [3l We carried out a search for the best-fit Keplerian orbital solution using a nonlinear, 
least-squares algorithm. In the fitting procedure, each velocity measurement is assigned a 
weight constructed from the quadrature sum of the internal measurement uncertainty and 
a stellar "jitter" term. The internal uncertainty is the weighted standard deviation of the 
mean velocity measured from each of the ~ 700 2-A chunks in each echelle spectrum. The 
stellar jitter term accounts for noise in excess of the internal uncertainties due to astrophys- 
i cal sources such as rotationa l modulation of stellar surface features and stellar pulsation 
( ISaar et al.lll998l : IWrightll2005l ). For Sun-like dwarfs, the stellar jitter can be accurately pre- 
dicted based on the chromospheric diagnostic logi?^^, given by the ratio of the CallH&K 
flux to the stellar UV continuum, along with luminosity and B — V color. However, logi?^^ 
is not calibrated for faint M dwarfs because stars with such cool atmospheres lack a well 
defined continuum due to molecular line blanketing. We therefore estimate the stellar jitter 
of the M dwarfs in our sample by assuming a constant value for all stars with B — V > 1.3. 
We selected the 133 M stars in our sample with 12 or more observations and measured the 



jitter using the prescription of IWrightl (120051 ). For HIRES observations made prior to the 
2004 CCD upgrade we derive a median jitter of 3.5 m s~^, and 2.0 m s~^ for post-upgrade 
observations. 

Our search for a best fitting Keplerian to the 18 observations of GJ317 reveals a min- 
imum in a/xF f*^^ orbital period near P = 690 days. However, the root-mean-squared 
(rms) residual to this best fitting, single-planet model is 16.8 m s~^, which is much higher 
than the expected scatter of ~ 6 m s~^ from the internal errors and jitter. The resulting 
reduced \/y?= 1.63 indicates the single-planet model is inadequate. 



X2 
; 

We attempted to improve on the single-Keplerian model by adding a variable linear 
trend. The existence of a trend in the velocities implies an additional, long-period orbital 
companion in the system. The dashed line in Figure [3] shows the best fitting Keplerian with 
a linear trend of 7.6±1 m s~^ yr~^. The inclusion of a trend reduces the rms from 16.8 m s~^ 
to 12.5 m s~^, and decreases a/xE^ from 1.63 to 1.42, after accounting for the extra free 
parameter. The Keplerian parameters of this fit are listed in Table [3l along with their 
estimated uncertainties. The parameter uncertainties were derived using a Monte Carlo 



method (e.g. iMarcy et al.l 120051 ). Our best-fit orbital solution yields a 1.897 year period, 
velocity semi-amplitude K = 71.0 m s~^, and eccentricity 0.193 ± 0.06. Using our adopted 
stellar mass M^, = 0.24 Mq, we calculate a minimum planet mass Mp sin z = 1.2 Mjup and 
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Fig. 3. — Our radial velocity time series for GJ317, based on observations from Keck observatory. The 
dashed line shows the best-fit Keplerian with an additional linear trend of 7.6 ± 1 m s^^year^^. The 
rms scatter of the residuals (bottom panel) is much larger than the predicted 6 m s^^ of scatter due to 
mearurement uncertainties and stellar jitter. 



semimajor axis a = 0.95 AU. 



3.1. The Two-Planet Model 

In light of the success of the Keplerian-plus-trend model, we performed a search for a 
double-Keplerian fit. The best fitting, two-planet model is shown in Figure HI with orbital 
periods of 673.4 and 2700 days, respectively. The minimum mass of the inner planet changes 
from 1.2 Mjup to 1.3 Mjup, compared the Keplerian-plus-trend model. The rms scatter 
drops from 12.5 m s~^ to 6.32 m s"^ with a resulting decrease in from 1.42 to 1.11. 

While the quality of the double-planet fit is encouraging, it is still possible that ran- 
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Fig. 4. — The best fitting two-planet model. The top two pannels show the single-planet fits with the 
other component removed. The addition of an outer planet with a period P w 2700 days decreases the rms 
scatter of the residuals from 12.5 m s~^ to 6.32 m s"^, and the reduced \/xi from 1.42 to 1.11. However, 
the duration and time-sampling our observations only marginally constrain the 11 free parameters of the 
two-planet fit. 



dom variability and jitter could conspire to produce a false curvature in our sparse series of 
measurements. If this is the case, then the second Keplerian and its associated 4 additional 
parameters are an unjustified embellishment on the single-Keplerian-plus-trend model. We 
tes ted this null hypoth esis by employing the scrambled-residuals false-alarm test described 
by I Wright et al.l (120071 ). We first subtract the single-Keplerian-plus-trend model from the 
data and adopt the residuals as random variates to construct a large number of mock ve- 
locity time series. Each mock time series is constructed by scrambling the residuals, with 
replacement, using a pseudo-random number generator. The scrambled residuals are then 
added back to original data to generate a mock data set. For each mock time series we 
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Fig. 5. — Empirical assessment of the false-alarm probability of the two-planet model. The observed 
velocities yield a double-planet orbital solution with \f~)(^= 1.11 (arrow). The histogram shows the distri- 
bution of \px^ obtained from singie-Keplerian fits to the mock velocity sets, with a variable linear trend. 
The mock data sets are generated by adding the scrambled residuals from a singie-Keplerian fit to the ob- 
served velocities. 11 of the 1000 mock data sets produced orbital fits with -^/xE^ lower than the original time 
series, indicating a FAP of 0.011 that random fluctuations are a viable explanation for the quality of the 
two-Keplerian fit to the observations. Additional observations are required to fully characterize the outer 
companion and reduce the FAP to below 0.01. 



perform a full search for a best-fit Keplerian-plus-trend orbital solution. If the curvature in 
the original data is an artifact of random noise, then many of the mock data sets should pro- 
duce lower a/xE^ than the two-planet fit to the original data. However, if a second, resolved 
Keplerian signal is present in the observations, then the mock data sets should produce fits 
that are much worse on average than the original velocities. 

The results of our test are shown in Figure O We generated 1000 mock data sets and 
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found that only 11 produced \fx^ equal to or lower than the original velocities, resulting in 
an false-alarm probability FAP= 0.011. The 1.1% FAP is promising, but we are nonetheless 
cautious about the validity of this two-planet model. Given the relatively low amplitude 
of the second signal and long period of the two candidate planets, our 18 measurements 
provide only 12 independent data points. The two-planet model is therefore only marginally 
constrained by our present data. Further, the duration of our observations only covers one 
orbit cycle of the putative outer companion. At this point we can only confirm the inner 
planet, and additional observations are required to verify the exact orbital properties of 
GJ317C. 



4. The Stellar Mass— Exoplanet Correlation 

4.1. Sample Selection 

The mass range encompassed by the California and Carnegie Planet Search target stars 
provides a large, uniform sample of targets that are ideal for evaluating the effect of stellar 
mass on the occurrence of Jovian planets. This analysis can be accomplished by simply 
measuring the fraction of stars with planets in three mass bins: the low-mass M and late 
K dwarfs with < 0.7 Mq, the Sun-like FGK stars with 0.7 < < 1.3 Mq, and the 
intermediate-mass subgiants with 1.3 < < 1.9 Mq. 

We were careful to select planets that could be reasonably detected around stars with 
a wide range of masses, and for surveys with different time baselines. Our first requirement 
is that stars have at least 8 observations spanning enough time to search for planets out to 

— 1/2 

a a distance of 2.5 AU. Using Kepler's third law, the orbital period scales as cx M," ' for 
a fixed semimajor axis. An orbital distance of 2.5 AU roughly corresponds to periods equal 
to the 3 year duration of the subgiants planet search, as well as the ~ 7 year duration of the 
NASA Keck M Dwarfs Survey. Next, we exclude stars sele cted as part of the metallicity- 



biased Next 2000 Stars (N2K) survey (IFischer et al.ll2005l ). Metallicity is an established 



tracer of giant planets, and we want to isolate the effects of stellar abundance from the 
effects of stellar mass. We also excluded stars that were added to the CCPS target list after 
a planet had already been announced by another group. 

We included only planet detections that induce a velocity semi-amplitude, K, large 
enough to be detected around stars in all three mass bins. For a planet of a given minimum 
mass and orbital period, the host star's velocity amplitude scales as K oc M^"^^^ . The larger 
masses of the subgiants the r efore limit the minimum planet mass that can be included in 



our analysis. iJohnson et al.l (120071 ) find that their Doppler measurements of subgiants yield 
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a typical precision of ~ 6 m s~^. We therefore include in our analysis only planets with 
minimum masses Mpsini > 0.8 Mjup. A planet with this minimum mass in a P = 3 year 
orbit would represent a 2-a detection around a 1.6 MQSubgiant, and would be much easier 
to detect around a lower-mass star. 

Our choice of such a large minimum planet mass is conservative since lower-mass planets 
would be readily detectable in shorter orbits. However, restricting our analysis to planets 
with masses greater than 0.8 Mjup also accomplishes a separate goa l. We wish to tes t the 
effects of stellar mass on the formation of giant planets. The models of llda fc Liru (j2005al ) and 
Laughlin et al.l (120041 ) demonstrate that the formation of planets like Jupiter is a threshold 
process limited by the growth rate of the rocky, embryonic cores. The core-growth rate is 
in turn limited by the supply of solid material at the midplane of the protoplanetary disk, 
which should scale with the mass of the central star. 



4.2. Jovian Planet Occurrence 



The results of our analysis are displayed in Figure [6] and are summarized here. For the 
low-mass bin, we selected stars with masses M^, < 0.7 M©. Of the 147 low-mass stars in the 
NASA M Dwarf Survey, 130 m eet our selection criteria. These stars have masses estimated 
using the iDelfosse et al.l (120001 ) K-band mass-luminosity calibration. We also calculated the 
masses of all stars with B — V > !.?> using the Padova stellar interior models (IGirardi et al. 



2OO2I ). This analysis revealed an additional 39 stars in the mass range 0.6 < < 0.7 M©. 
Of these 169 total stars, only 3 are known to harbor at least one giant planet: G J 876, G J 849 
and GJ317. The resulting planet occurrence rate for these low-mass stars is 1.8 ± 1.0%. 

For the Solar-mass bin we selected stars from the CCPS Lick, Keck and Anglo- Australian 
Observatory surveys with masses in the range 0.7 < < 1.3 M©. A total of 803 stars met 
our stellar mass and planet-detectability criteria. Of these stars, 34 harbor at least one giant 
planet, resulting in a planet occurrence rate of 4.2 ± 0.7%. 



For the high-mass bin, we used the intermediate-mass subgiants in the iJohnson et al. 



( l2006al ) sample, of which 68 me et our selection crite ria and have masses in the range 1.3 < 
< 1.9 Mq according to the iGirardi et al.l (120021 ) stellar interior models. We found an 
additional 33 intermediate-mass subgiants in the regular CCPS Keck and Lick samples, 
for a total of 101 stars. There are 6 published planetary companions among these stars. 
An additional 3 companions will be announced in a future publication pending follow-up 
observations at Lick Observatory in the upcoming observing seasons (Johnson et al. 2007b, 
in preparation). While these additional planet candidates lack sufficient phase coverage 



-14- 



to publish at this time, their signals currently have false-alarm probabilities FAP< O.Olo. 
Including these strong candidates results in 9 giant planets among 101 intermediate-mass 
subgiants, for an occurrence rate of 8.9 ± 2.9%. 



4.3. The Effects of Metallicity 



The occurr ence of detecta b le planets has previously been shown to co rrelate with stel- 
lar metallicity (iGonzalea 119971 : ISantos et al.ll2004l : iFischer fc Valentil l2005l ) . It is therefore 
important to test whether the correlation between stellar mass and planet occurrence shown 
in Figure E] is real, or merely an effect of an underlying metallicity bias. 

We searched for such a bias in our sample by measuring the median metallicity in each 
mass bin. Since the Solar-mass stars and intermediate-mass subgiants were analyzed using 
the same LTE spectral synthesis code (SME), our results from these two mass regimes can be 
directly compared. For the Solar-mass bin we find a median metallicity of [Fe/H]= +0.042. 
The median metallicity of the subgiants is slightly higher at [Fe/H]= +0.075. In order to 
compare the expected occurrence rate in each of these mass bins based on metallicity alone, 
we ena ploy the relationship between metallicity and planet probability from lFischer &: Valenti 
fl2005h . 



P(planet) =0.03 X 10'[^^/«]. (1) 

Given the median metallicity of the Solar-mass stars. Equation [T] predicts a 3.6% prob- 
ability of finding a detectable planet. This is consistent within errors with the 4.2 ± 0.7% 
probability that we measure. For the subgiants, the expected probabihty is 4.2%, which is a 
factor of 2 lower than the 8.9 ± 2.9% we measure. This enhanced planet occurrence around 
our high-mass stars is significant at the 1.6-0" level. 

The correlation between stellar mass and exoplanets in Figure [6] is most apparent when 
comparing the bins at either mass extreme. The high-mass stars in our sample appear to 
be planet-enriched by a factor 5 compared to the low-mass M dwarfs, which represents a 
2.3-cr result. However, it is currently very difficult to derive spectroscopic abundances for M 
dwarf^. Could it be that there is a systematic metallicity bias among our low-mass stars? 



^We have analyzed the velocities for our entire sample of M dwarfs and find no candidate signals due to 
Jovian-mass planets with FAP< 0.01. Similarly, excluding the metal-rich N2K targets, none of the FGK 
stars in the CCPS samples has an unpublished candidate with Mp sini > 0.8 Mjup within 2.5 AU. 

^See, however. iBean et al.l (|2006l ) for spectroscopically derived metallicities of 3 M dwarf planet host stars. 
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Stellar Mass [M©] 

Fig. 6. — Histogram illustrating the rising percentage of stars with detectable planets as a 
function of stellar mass. The stars selected for each mass bin have 8 or more observations 
providing detectibility of planets with masses Mp sin i > 0.8 Mjup out to a = 2.5 AU. After 
correcting the measured percentages in each mass bin for the effects of stellar metallicity, the 
rising trend is slightly diminished (filled circles, offset to the right by +0.05 M0 for clarity). 
However, the high-mass bin is uncorrected for the decreased sensitivity of Doppler mea- 
surements of higher-mass subgiants compared to lower-mass stars (cf § 14.31) . Our measured 
occurrence rate for high-mass stars therefore represents a lower limit unlike the Solar-mass 
and low-mass bins. The error bars on each bin are from Poisson statistics and the numbers 
above each bin compare the number of stars with planets N^osts to the total number of 
stars in each bin N stars- 
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A c ommon argument fo r such a metallicity bias among M dwarfs is related to their 



age (e.g. iBonfils et aLll2005bl ). According to this argument, M dwarfs have hfetimes longer 
than the age of the Galaxy and are therefore more likely than higher-mass stars to have 
formed in the distant past when the Galaxy was metal poor. This argument is based on the 
assumption that there is a well-defined age-metallicity relationship in the Galaxy. However, 
studies of Galactic chemical evolution have shown that no well-defined correlation between 
age and metallicity exists in the thin Galactic disk. Analyses of the local stellar population 
show that stars have a wide range of abu ndances at ages from 1 to 10 Gyr, with a 0.2 



dex scatter in [Fe /H] over this age range (lEdvardsson et al.l Il993l : iNordstrom et al.l 12004 



Takeda et al.ll2007l ). Similarly, open clusters with comparable ages also exhibit a wide range 



of me tallicities, with 3 to 6 Gyr clusters differing by as much as 0.82 dex (e.g. iChen et al. 
20031 ) . Further, late G- and all K-type stars also have lifetimes comparable to the age of 
the Galaxy. Therefore, any metallicity bias present in our sample of M dwarfs should also 
be present in a volume-limited sample of G and K dwarfs. We therefore find no compelling 
reason to suspect that the abundances of M dwarfs should deviate significantly from the G 
and K stars in the Solar neighborhood. 

Our sample of M dwarfs within 15 pc of the Sun should trace the metal abundances 
of the other stars in th e immediate Solar neighborhood. The volume-limited sample of 
Fischer &: Valentil (120051 ) has a mean metallicity of [Fe/H]= —0.05. This is roughly the 
abundance we expect for our sample M dwarfs. We stress that, contrary to the argument 
above, our sample of M dwarfs is not metal-poor compared to other nearby stars. Rather 
it is the Solar-mass stars and subgiants among the planet search targets that are slightly 
metal rich compared to a volume-limited sample. This effect has been previously noted by 
Fischer fc Valentil (120051 ) and is related to the fact that in a given range of B — V colors 
metal-rich stars are brighter a nd more likely to be included in the samples of most Doppler- 
based planet searches (see also lMarcy et al.ll2005l . for a more detailed discussion.). 



Based on our estimate for the mean metallicity of our sample of M dwarfs ([Fe/H] = 
—0.05), Equation [1] predicts a 2.4% planet occurrence rate for the stars in our low-mass bin. 
This prediction is slightly higher than our measured occurrence rate of 1.8 ± 1.0%, but the 
two values agree within errors. 

We can use the probabilities predicted by Equation [T]to correct our measured occurrence 
rates for the effects of stellar metallicity. Each corrected planet fraction, //, is constructed 
as 



fl = (PQ/Pi) X fi, 



(2) 
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where is the original planet fraction, Pq = 0.03 is the expected probability for solar 
metallicity, and Pi is the predicted probability from Equation [TJ Similarly, the corrected 
uncertainty in each bin is 



X (Ti 



(3) 



The filled circles in Figure [H] show the corrected values and uncertainties of the planet oc- 
currence rate in each mass bin. After correcting for metallicity, the contrast between the 
low-mass and high-mass bins is diminished. However, evidence for an increasing trend 
remains: the A-type stars have a factor of 2.5 higher planet occurrence rate than the M 
dwarfs. 

It is important to note that the corrected occurrence rates in Figure [6] do not account 
for the decreased sensitivity of our Doppler survey of subgiants compared to the surveys of 
lower-mass stars. Our requirement of 8 or more observations was exceeded in most cases 
by the stars in the Solar-mass and low-mass bins. The Solar-mass stars had on average 
21 observations each, and the M dwarfs averaged 16 observations. On the other hand, the 
subgiants had an average of only 10 observa tions each. Also, whereas our typ i cal Do ppler 
precision is <2 m s~^ for Solar-mass stars (IMarcy et al.l l2005l : iJohnson et al.l l2006bl ) and 
3-4 m s~^ for M dwa rfs (Figure ED, the su bgiants typically have a precision of 4-7 m s~^ due 
to their larger jitter (1 Johnson et al.ll2007l ). The sensitivity of our measurements of subgiants 
is further diminished since the amplitude of a star's reflex motion caused by a planet with a 

2/3 

given mass and orbital period scales as cx M* . We therefore conclude that the measured 
planet occurrence rate for the stars in our high-mass bin represents a lower limit rather than 
an absolute measurement of the true fraction of subgiants with planets (as indicated by the 
upward-facing arrow in Figure E]), even after correcting for metallicity. 



5. Summary and Discussion 



We present the detection of a Jupiter-mass planet orbiting the M3.5 dwarf GJ 317. This 
detection marks the sixth M dwarf known to harbor at least one Doppler-detected planet. 



toget her with G J876 flMarcv et al.lll998l : lOelfosse et al 



1998 



Marcv et al 



200 



2005h. GJ436 JSutler et al.l l2004h. GJ581 JBonfils et al.l l2005bl: ludrv et al 



Rivera et al. 



2003), GJ849 



flButler et al.ll2006ar ) and GJ674 flBonfils et all 120071 ). GJ317 is the only third M dwarf out 
of roughly 300 surveyed known to harbor a Jovian planet. In addition to our firm detection 
of a Jovian planet with P = 692.9 days, we also detect evidence of a possible second Jovian 
planet in the system near P ^ 2700 days. However, additional monitoring is required to 
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fully characterize the orbit of the second planet. 



Multi-planet systems appear to be relatively common among M dwarfs compared to 
Sun-like stars. All M stars with one Jovian planet show evidence of a second companion. 
GJ876 has a pair of Jupiter-mass planets in a 2:1 mean motion resonance, along with an 



inner "super Earth" (IMarcv et al.ll200ll:lRivera et al.ll2005r) . GJ849 has a long-period Jovian 



planet with a linear trend (IButler et al.l l2006al ). Of the 3 M dwarfs with Neptune-mass 
planets, two have m ultiple planets or evidence of an additional compar iion: GJ581 harbor s 
3 low-mass planets (lUdry et al.l 120071 ). and GJ436 has a linear trend (IManess et al.l 120071 ). 
Only GJ674 appears to be in a single-planet system. From the first 6 planet detections 
around low-mass stars, it appears as though M dwarfs have an 80% occur rence rate of 



multi -planet systems, compared to the 30% rate measured for FGK stars (IWright et al. 
2007h . 



The high frequency of multi-component detections around M dwarfs may be in part 
due to the increased detectability of planets around low-mass stars. This is likely the case 
for the triple Neptune system around GJ581. The "c" component in that system would 
have induced a velocity amplitude of only 1.4 m s~^ if it orbited a 1 Mq star. However, the 
enhanced detectability of planets around M stars shouldn't matter as much for the systems 
containing Jovian planets since all of the giant planets orbiting M stars would be detectable 
around Solar-mass stars. 

We studied the relationship between stellar mass and the occurrence rate of giant plan- 
ets by combining our sample of M dwarfs with our samples of Solar-mass FGK stars and 
intermediate-mass subgiants. By measuring the fraction of stars with planets in three stel- 
lar mass bins we find that the frequency of planets with > 0.8 Mjup within a < 2.5 AU 
increases with stellar mass (Figure E]). The retired A-type stars in our sample are nearly 5 
times more likely than M dwarfs to harbor a giant planet. This important result establishes 
stellar mass as an additional sign post for exoplanets, along with metallicity. Just as metal- 
licity informs the target selection of searches for short-period planets (e.g. iFischer et al. 



20051 ). stellar mass will be an important factor in the target selection of future high-contrast 
direct imaging surveys. While the lower luminosities of M dwarfs provide favorable contrast 
ratios that facilitate the detection of thermal emission from young giant planets, our results 
show that A-type stars are far more likely to harbor such planets. 

In order to understand the role of stellar mass on planet formation, it is important to 
disentangle the known effects of stellar metallicity from our stellar sample. After correcting 
for metallicity in each mass bin, the slope of the trend in Figure E] is slightly diminished. 
However, the factor of 2.5 increase in planet occurrence around A stars compared to M 
dwarfs remains significant, especially considering that the high-mass bin is uncorrected for 
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the relatively lower detection sensitivity of the subgiants planet search. Our results therefore 
confirm the prediction o f the core accretion naodel that higher raass stars should forna Jovia n 
planets more efficiently (iLaughlin et al.ll2004j : llda &: Lirul2005bl : Kennedy &: KenyonI 120071 ) . 



Our discovery confirms and expands upon the recent results of iLovis fc Mayorl (120071 ). 
who are searching planets around intermediate-mass K giants. By focusing on stars in open 
clusters, they are able to control for the age and metallicity of their sample and determine 
accurate stellar masses. Because of the larger stellar radii and lower precision (jitter > 
15 m s~^) of K giants, they focused on planets with Mp sini > 5 Mjup and 0.5 > a > 2.5 AU. 
Lovis &: Mayorl (120071 ) find that the occurrence rate of "super-Jupiters" and brown dwarfs 
within this mass and semimajor axis range increases with stellar mass, rising from 0% for M 
dwarfs to 2.5% for evolved A stars. 



The high-mass sample of stars analyzed by iLovis &: Mayorl (120071 ) contained primar- 
ily stars with masses M^, > 2.0 M©. In order to avoid confusion with horizontal-branch 
stars ("clump giants") we have restricted our subgiants planet search primaril y to stars with 
absolu te magnitudes My > 1.8, corresponding to masses < 2.1 M© (I Johnson et al. 
2006br) . Our planet search is therefore complement ary to the survey of iLovis fc Mavor 



J2007h and other p l anet searches around K giants (e.g. ISetiawan et al.ll2003l : ISato et al 



20031 : 



Hatzes et al.ll2005l : iReffert et al.ll2006l : iNiedzielski et al.ll2007n . The K giants provide infor- 
mation about massive planets around stars with M* > 2.1 M©, and the relatively stable 
atmospheres of subgiants (jitter = 5 m s~^) allow us to detect planets beyond 1 AU with 
minimum masses down to Mp sinz ^0.5 Mjup. 



Our preliminary results and those of iLovis fc Mayoii (120071 ) reveal a positive correlation 
between the mass of stars and the likelihood that they harbor giant planets. However, 
testing the exact shape of the relationship predicted by the simulations of lKennedy fc Kenyon 
(120071 ) — a rising trend of planet occurrence up to a peak near 2.5 Mq and a decreasing 
trend toward higher masses — will require a much larger sample of stars. We have expanded 
our search for planets around intermediate-mass stars by adding 300 additional former A 
and F stars to our Lick and Keck samples. The results from our expanded planet search 
should reduce the error bar on the high-mass bin in Figure [HI by a factor of 2. If the 9% 
occurrence rate for > 1.3 M© is confirmed, then our expanded planet search will result 
in the detection of 20-30 new planets orbiting some of the most massive planet host stars 
in the Solar neighborhood. In addition to verifying the preliminary results presented here, 
these planet detections will also allow us to study the effects of stellar mass on other planet 
characteristics such as eccentricity, semimajor axis and minimum mass. 
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Table 1. Radial Velocities for GJ 317 



JJJ 


rlV 


Uncertainty 


Oil r\r\r\r\ 
-Z [ iUUIJU 


( ill S j 


(ill S j 


11550.993 


13.89 


3.79 


11552.990 


30.27 


4.99 


11582.891 


45.77 


4.16 


11883.101 


-23.35 


3.74 


11973.795 


-65.89 


6.10 


12243.073 


0.00 


6.52 


12362.949 


95.56 


6.07 


12601.045 


-31.64 


5.09 


12989.125 


101.25 


4.79 


13369.016 


-29.65 


3.40 


13753.983 


130.31 


3.68 


14084.001 


-6.55 


4.64 


14086.141 


-21.73 


4.55 


14130.082 


-15.26 


4.79 


14131.014 


-15.47 


4.20 


14138.932 


-8.98 


2.82 


14216.733 


4.15 


3.86 


14255.745 


28.54 


2.25 
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Table 2. Stellar Properties for GJ 317 



Parameter 


Value 


V 


12.0 


B-V 


1.52 


Mv 


12.2 


K 


7.016 


J-K 


0.915 


H-K 


0.311 


d (pc) 


9.17 ±1.7 


(Mo) 


0.24±0.04 


[Fe/H] 


-0.23±0.2 


Table 3. Orbital Solution for GJ 317 b 


Parameter 


Value 


P (days) 


692.9±4 


P (years) 


1.897±0.01 


K (ms-^) 


71.0±7 


e 


0.193±0.06 


Tp (Julian Date) 


2451639±30 


Linear Trend (m s"^ 


year-i) 7.6±1 


Lo (degrees) 


344±10 


Msini (Mj^p) 


1.2 


a (AU) 


0.95 


Nobs 


18 


RMS (ms-^) 


12.5 




1.42 



